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THE VIBRATION CHARACTERISTICS OF A COUPLED HELICOPTER 


ROTOR-FUSELAGE BY A FINITE ELEMENT ANALYSIS 

Michael J . Rutkowski 

Ames Research Center 
and 

Aeromechanics Laboratory 
AVRADCOM Research and Technology Laboratories 


SUMMARY 


The dynamic coupling between the rotor system and the fuselage of a simplified 
helicopter model in hover is analytically investigated. Mass, aerodynamic damping, 
and elastic and centrifugal stiffness matrices are presented for the analytical 
model; the model is based on a beam finite element, with polynomial mass and stiff- 
ness distributions for both the rotor and fuselage representations. For this 
analytical model, only symmetric fuselage and collective blade degrees of freedom 
are treated. Real and complex eigenanaly ses are carried out to obtain coupled 
rotor-fuselage natural modes and frequencies as a function of rotor speed. Vibra- 
tion response results are obtained for the coupled system subjected to a radially 
uniform, harmonic blade loading. The coupled response results are compared with 
response results from an uncoupled analysis in which hub loads for an isolated 
rotor system subjected to the same sinusoidal blade loading as the coupled system 
are applied to a free-free fuselage. It is shown that although the fuselage 
response as a function of forcing frequency and rotor speed is similar in the two 
cases, the responses resulting from the approximate, uncoupled analysis are signifi- 
cantly greater than those resulting from the coupled analysis. Thus, it is necessary 
to carry out a coupled rotor-fuselage analysis in order to accurately predict the 
fuselage vibration response. 


INTRODUCTION 


Vibration is a serious problem in military and civilian helicopters. It causes 
structural fatigue and failure, increases maintenance and operational costs, impairs 
performance, and has important human factor implications. Surprisingly, however, 
vibration reduction efforts during the original design of most helicopters are not 
significant. Too often it remains for later and costly flight-test programs to 
identify vibration problems and to suggest fixes — for example, the addition of vibra- 
tion absorbers or, if necessary, even structural modifications. In addition to the 
costs they impose, these ad hoc fixes usually carry with them a weight penalty that 
reduces either payload or range or both. 

There are many sources of vibration excitation in a helicopter — for example, 
engine and gearbox vibratory forces, tail-rotor excitation, rotor mass imbalance, and 
rotor wake impingement — but main— rotor excitation is usually the most significant. 
This excitation results from a combination of phenomena such as hover or forward 
flight aerodynamics, rotor inflow, hub motion, blade dynamics, and coupled rotor- 
fuselage dynamics. Unfortunately, the complete helicopter vibration problem is 



virtually impossible to predict, as yet, because of its structural, as well as aero- 
dynami c , complex! ty . 

To make the helicopter vibration problem more tractable, simplifying assumptions 
must be made. Typically, the rotor system and the fuselage are analyzed separately 
and then attempts are made to account for rotor-fuselage coupling. First, hub forces 
and moments are obtained for an isolated elastic rotor rigidly supported at the hub 
and subjected to aerodynamic loads. The hub loads are then applied to a flexible 
fuselage, such as a NASTRAN finite element model, and responses are then calculated 
at different points on the fuselage. This procedure is a questionable one, however, 
because it does not account for rotor-fuselage coupling. In reference 1, it is noted 
that even introducing an equivalent rotor mass, although it gives a somewhat better 
approximation, still does not adequately represent the coupled rotor-fuselage sys- 
tem. Thus, there appears to be no viable alternative to carrying out a coupled 
rotor-fuselage analysis when one is investigating fuselage response to rotor 
excitation. 

Recently, several coupled rotor-fuselage analyses were carried out using simple 
mathematical models in order to gain insight into the physics of the coupled vibra- 
tion problem (refs. 1-4). In reference 1, a coupled rotor-airframe analysis based 
on an impedance matching method is described. This analysis is then applied to a 
simple, single-rotor example and to predicting vibration in the Model 347 helicopter. 
The role of rotor impedance in rotorcraft vibration is assessed in reference 2 by 
Hohenemser ' and Yin. In particular, they show that the usual approximation whereby 
dynamic forces and moments for a rigidly supported isolated rotor are applied to a 
flexible airframe can lead to large errors in vibration predictions. In their 
improved method, Hohenemser and Yin use the rotor impedance to correct the dynamic 
rotor force and moment inputs to the airframe. This work was extended by Hsu and 
Peters (ref. 3) who performed a coupled rotor-airframe analysis using harmonic bal- 
ance by matching fuselage impedance with rotor impedances (flapping only) calculated 
for a wide range of rotor parameters. The results of these analyses show that hub 
motions have their greatest effect on the hub loads of relatively stiff rotors. A 
harmonic balance solution was also used by Kunz (ref. 4) in a fully coupled rotor- 
craft vibration analysis which showed the significant effects of rotor-body coupling 
upon hub pitching moment. 

The objective of the present work is to extend the previous analyses (refs. 1-4) 
by carrying out a qualitative examination of both the basic rotor-fuselage coupling 
mechanism and the effect of this dynamic coupling on fuselage response. This report 
presents the results of an investigation into the dynamic coupling between the rotor 
system and the fuselage of a simplified beam-model helicopter. The rotor-fuselage 
model used in this analysis is based on a two-degree-of-f reedom (transverse deflec- 
tion and rotation) beam finite element with polynomial mass and stiffness properties. 
Both real and complex eigenanalyses are carried out to obtain the symmetric fuselage 
modes and the blade collective modes as a function of rotor speed. Mode shapes, as 
well as modal frequencies, are obtained in both the real and complex analyses. In 
the complex eigenanalyses, damping for the rotor in hover is introduced using quasi- 
steady strip-theory aerodynamics. In addition to these eigenanalyses, this report 
also describes forced response vibration results for the case of a radially uniform, 
harmonic load applied to the rotor. Results from coupled, as well as uncoupled, 
rotor-fuselage response analyses are presented and compared to demonstrate the impor- 
tance of dynamic coupling. These results are presented in several different formats 
including two- and three-dimensional plots and contour plots. 
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Because of the relative simplicity of the model used in these analyses, the 
present work considers only symmetric motion of the rotor and fuselage. Nonetheless, 
there is still a great deal that can be learned from the results of these analyses. 
The model^s simplicity allows one to easily generate responses at any node point for 
a large number of forcing frequencies at many rotor speeds. 


MATHEMATICAL MODEL 


The basic structural element used to model the rotor and the fuselage in this 
analysis is the beam-bending finite element, as described in reference 5. For sim- 
plicity, this element (fig. 1), which is based on a cubic polynomial displacement 
function, includes transverse displacements and rotations in one plane only. The 
element mass, damping, and stiffness matrices used in this analysis are derived in 
appendix A and presented in appendix B. The usual uniform elastic stiffness and con- 
sistent mass matrices have been extended to allow for up to a cubic polynomial varia- 
tion in both the element mass and stiffness properties. In addition, the rotor-blade 
elements also include centrifugal stiffness matrices and aerodynamic damping matrices 
based on quasi-steady strip theory. Both of these blade element matrices vary with 
rotor speed. Structural damping is included in both the rotor and the fuselage 
elements . 

Because of the relative simplicity of the beam finite element used to model the 
rotor and fuselage, the present investigation considers only the rotor-induced ver- 
tical shear loads that are transmitted to the fuselage. The longitudinal and lateral 
forces, which cause pitching and rolling in an actual helicopter, are not included 
here. Therefore, this coupled analysis is restricted to symmetric loadings and model 
configurations so that only the collective rotor modes and symmetric fuselage bending 
modes are excited. 

The equations of motion for the real and complex eigenanalyses and the forced- 
response analyses are given in appendix C, In the forced-response analyses, in addi- 
tion to the aerodynamic loading, the rotor system is subjected to a distributed har- 
monic load. Although in the present analysis a radially uniform load was considered, 
a load distribution that varies linearly or quadratically with blade span, as well as 
one or more point-loads applied to the rotor or the fuselage or both, could be con- 
sidered. The present investigation can readily be extended to a study of the sensi- 
tivity of coupled rotor-fuselage vibrations to small variations in geometric, inertial, 
and stiffness parameters. In addition, natural antiresonances, as well as the effect 
of introducing an isolator between the rotor and the fuselage, could also be 
investigated. 

The rotating beam mathematical model used in this analysis has been validated 
for an isolated rotor blade in reference 6. Real eigenvalue and eigenvector results 
obtained using the beam element and a cubic polynomial displacement function were 
compared with the nearly exact results obtained using a higher-order finite element 
for several rotating beam examples. 


RESULTS AND DISCUSSION 


The beam model representation of the basic configuration used in this coupled 
rotor-fuselage analysis is shown in figure 2. This model consists of a hingeless 
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rotor — made up of two 5-element rotating beams cantilevered in bending to a lumped 
mass hub connected by a linear spring to the midpoint of a 10-element fuselage. 

A 5-element representation of the blades was chosen because it has been shown to 
give frequencies within 1% of the exact values for the first three rotating modes 
(ref. 6). Figure 2 also shows the location of the rotor and fuselage node points for 
the present configuration. The blades are each 7.62 m (25.0 ft) long with uniform 
mass and stiffness distributions of 14.17 kg/m (0.30 slugs/ft) and 82,600 N«m^ 

(200,000 lb*ft^), respectively, and the concentrated hub has a mass of 233.5 kg 
(16.0 slugs). In the present configuration the connection between the rotor and the 
fuselage is assumed to be essentially rigid, and, therefore, the connecting spring 
has a relatively high spring constant, k, of 1.46x10^^ N/m (1.00x10^^ Ib/ft). The 
flexible fuselage is 12.2 m (40.0 ft) long with the nonuniform, but symmetric, mass 
and stiffness distributions shown in figures 3 and 4, respectively. These mass, 
stiffness, and geometric properties were chosen to demonstrate some of the dynamic 
coupling characteristics of coupled rotor-fuselage systems and still to give frequency 
placements that are representative of typical helicopters having a nominal operating 
rotor speed of = 30 rad/sec. In addition, a reasonably typical rotor-plus-hub mass 
to fuselage mass ratio of 0.11 was chosen. Finally, in the response analyses, in 
addition to the aerodynamic damping in the blades, structural damping of 6% has been 
assumed throughout the rotor-fuselage in order to reduce the magnitudes of the reso- 
nance peaks in the response plots. 


Eigenanalyses 

Natural frequencies— The uncoupled and coupled modal frequencies as a function 
of rotor speed are presented in figure 5(a) for the rotor and fuselage without struc- 
tural or aerodynamic damping. Also shown in the figure are the (dotted) nP lines 
(n = 1,2, 3, 4) where the modal frequency is equal to the rotor speed or one of its 
harmonics. The boundary conditions for the uncoupled systems were taken to be free- 
free for the fuselage and hub-fixed for the rotor, since these boundary conditions 
most closely approximate the actual boundary conditions of each component in the 
coupled system for the typical rotor-mass/fuselage-mass ratio considered. For the 
uncoupled free-free fuselage, figure 5(a) shows the two symmetric bending frequencies 
(Fi and F^) to be independent of rotor speed. The uncoupled, collective flap rotor 
modes (B^ , B2 , and B3), on the other hand, increases in frequency as rotor speed is 
increased and, consequently, as centrifugal stiffening is increased. Furthermore, 
when the fuselage and rotor are uncoupled, the curves of rotor flap and fuselage bend- 
ing frequencies are seen to cross as the flap stiffness increases. Similar frequen- 
cies are obtained for the coupled rotor-fuselage system except where the rotor and 
fuselage couple. When the fuselage and rotor are elastically coupled, the frequencies 
are close to their uncoupled values only when the frequencies are widely separated. 
When a region where the uncoupled frequencies cross is approached, the frequency 
curves for the coupled modes begin to diverge and approach the frequency curves of 
the companion modes. In the present paper this region of modal coupling is referred 
to as a region of modal coalescence. Here, for the in vacuo case, figure 5(a) shows 
that the uncoupled frequencies are poor approximations to the actual coupled frequen- 
cies, especially near where the second fuselage-bending mode couples with the second 
rotor flap mode. 

The effect of rotor aerod 3 mamic damping on the uncoupled rotor and the coupled 
rotor-fuselage system frequencies is shown in figure 5(b). Except for the reduction 
of the first rotor-f lap-mode frequency from above IP to below IP, rotor aerodynamic 
damping does not have an appreciable effect on the uncoupled rotor modes. For the 
coupled rotor-fuselage system, however, the addition of aerodynamic damping results 
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in a significant change in the shape of the curves in the region of the rotor- 
fuselage mode frequency crossings. Figure 5(c) shows that in the presence of aero- 
dynamics, the effect of rotor-fuselage coupling is less than it is in the in vaouo 
case. Therefore, with rotor aerodynamic damping included, the uncoupled rotor and 
fuselage frequencies are a somewhat better approximation to the coupled system 
frequencies. 

Mode shapes— Mode-shape plots are presented to better illustrate the phenomenon 
of modal coupling and, in particular, to show what happens when rotor and fuselage 
modes couple. In addition, mode-shape plots are useful in showing the often over- 
looked effects of aerodynamic damping and out-of-phase modal content on coupling. 

In any discussion of mode shapes one inevitably faces the problem of how to name 
or label different modes that couple and "cross” on the f requency-versus-rotor-speed 
plot. Some liberties must be taken at times, since at modal crossings individual 
mode identities may become lost or ambiguous. For example, modes could be labeled to 
represent the order of value of mode natural frequencies, that is, mode 1, mode 2, 
etc. For a given "physical mode," however, such a mode number changes at frequency 
coalescences; for example, the "first rotor-flap mode" starts off as mode 1 at low 
rotor speed and switches to mode 2 at high rotor speed. In this report, modes will 
be referred to by their dominant physical characteristic. 

The eigenvectors for the first two rotor collective flap modes are shown in fig- 
ures 6(a) and 6(b), and the first two fuselage symmetric bending modes are shown in 
figures 6(c) and 6(d). The mode shapes for these four modes are presented for each 
of three conditions; (1) either an uncoupled hub-fixed rotor in vaouo or an uncoupled 
free-free fuselage; (2) coupled rotor-fuselage in vacuo \ and (3) coupled rotor- 
fuselage with rotor aerodynamic damping included. The mode shapes for each case, 
except the two uncoupled fuselage modes, which are independent of rotor speed, are 
presented for a range of rotor speeds. For the coupled rotor-fuselage system with 
aerodynamic damping included, both the real and imaginary parts of the complex mode 
shapes are plotted for rotor speeds from 0 to 60 rad/sec. In most cases the imaginary 
parts of the mode shapes, which come entirely from aerodynamic and structural damping, 
are small; as a result, the mode shapes are essentially real and differ only slightly 
from the mode shapes that were obtained from a real eigenanalysis with damping 
neglected. 

Consider first the modes of the isolated rotor blade that are shown in fig- 
ures 6(a) and 6(b). Of these two uncoupled modes, which change somewhat with rotor 
speed, the change in the first rotor mode shape as rotor speed increases is most 
apparent (fig. 6(a)). Initially, at zero rotor speed, the mode shape of each blade 
is that of a uniform cantilever beam. Then, as rotor speed increases, the rotor- 
blade stiffness increases because of increased centrifugal forces, and the blade 
mode shape approaches a straight line. Similarly, the second uncoupled rotor mode 
also tends to flatten out as rotor speed is increased (fig. 6(b)). 

The first two symmetric bending modes for the uncoupled free-free fuselage are 
shown in figures 6(c) and 6(d). These modes do not change shape as rotor speed 
increases, since the uncoupled fuselage modes are independent of rotor speed. 

Now we will trace the changes in these four basic modes (first rotor-flap, second 
rotor-flap, first fuselage-bending, and second fuselage-bending) as they are mani- 
fested in the coupled system. Usually their identity is evident, but when modes 
coalesce their separate identities are often lost. The effect of rotor-fuselage 
coupling on the rotor-flap mode shapes is small. For the coupled rotor-fuselage 
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configuration in figures 6(a) and 6(b), both with and without aerodynamic damping 
included, the rotor deflections are nearly the same as for the uncoupled rotor-blade 
modes. At the same time, the fuselage deflections for these rotor modes are very 
small, except at rotor speeds at which there is a coalescence of a rotor mode with a 
fuselage mode. The rotor modes do not exhibit much fuselage motion, because the mass 
of the fuselage plus the concentrated hub mass is much larger than the mass of the 
rotor blades. (In this case the mass ratio is 20:1.) In addition, the hub motion 
is usually negligible; this justifies the conventional uncoupled approximation to the 
coupled rotor-fuselage system whereby hub forces and moments, which are obtained for 
an isolated elastic rotor rigidly supported at the hub, are applied to a flexible 
fuselage. In contrast, the fuselage modes (figs. 6(c) and 6(d)) generally produce a 
large hub motion which in turn “excites” large rotor deflections in fuselage modes. 

Of particular interest is the behavior of the coupled rotor-fuselage configura- 
tion near = 30 rad/sec, where the second fuselage bending mode is highly coupled 
with the second rotor-flapping mode. The in vacuo mode shapes at = 30 rad/sec 
for these two modes (figs. 6(b) and 6(d)) look very similar except that the fuselage 
deflections are 180® out of phase. Both of these real modes contain fuselage-bending 
and rotor-flapping deflections. The mode shapes in air for these two modes, on the 
other hand, are very different from each other. This difference shows up in the com- 
plex eigenvectors. In figure 6(b), for second rotor-flap mode in air, the real part 
of the modal deflections contains only second rotor-flapping, and the imaginary mode 
shape contains a very small amount of both rotor and fuselage deflections. In fig- 
ure 6(d), for the second symmetric fuselage-bending mode in air, however, the real 
part of the complex mode shape is predominantly second fuselage-bending, and the 
imaginary mode shape shows a very significant amount of second rotor-flapping. 

To better focus on the details in the = 30 rad/sec region of high rotor- 
fuselage coupling previously seen in the frequency plot of figure 5(c), mode shapes 
for rotor speeds of 24-36 rad/sec are presented in figures 7(a) and 7(b). These fig- 
ures show the changes that occur in the coupled rotor-fuselage second rotor- and 
second fuselage-mode shapes, respectively, in the vicinity of their coalescence. In 
figure 7(a), the fuselage deflections for the second rotor mode in vacuo are seen to 
first increase as the rotor speed is increased to 30 rad/sec and then to decrease 
with opposite phase as the rotor speed is increased further. In air, both the real 
fuselage deflections and all of the imaginary deflections are small. This implies 
that for the rotor-f lapping and fuselage-bending modes considered here, rotor aero- 
dynamic damping suppresses the effects of rotor-body dynamic coupling. In figure 7(b), 
however, for the second fuselage mode in vacuo^ the rotor deflections first increase 
and the fuselage deflections decrease as rotor speed is increased to approximately 
30 rad/sec. Then the rotor deflections change phase and decrease with further 
increase in rotor speed while the fuselage deflections increase. Again, the in-air 
deflections are considerably different. They show a gradual change in the real part 
of the rotor mode shape from second mode to first mode of opposite phase, while the 
imaginary part of the rotor deflections first increases and then decreases. A similar 
type of behavior occurs for the coalescence of the first rotor-blade flap mode with 
the first fuselage-bending mode. The in vacuo and in- air mode shapes for these two 
modes are shown in figures 7(c) and 7(d). Again, as in the previous instance, the 
in vacuo and in-air shapes are considerably different near the coalescence. 

As has just been seen, in following a mode for increasing rotor speed there can 
be a significant difference in behavior between the in vacuo and the in-air cases, 
especially near a coalescence. In figure 8 are shown the uncoupled, coupled in vacuo y 
and coupled in-air eigenvalues as a function of rotor speed for the two modes pre- 
sented in figures 7(a) and 7(b). The second fuselage-mode eigenvalues for the coupled 
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in-air case are seen to fall roughly midway between those for the free-free fuselage 
alone and for a free-free fuselage with an effective rotor mass that is an approxima- 
tion sometimes used to account for rotor-airframe interaction (ref. 1). Similarly, 
the second rotor-mode eigenvalues for the coupled in-air case roughly agree with those 
for an isolated rotor both with and without aerodynamic damping included. The two 
in vaouo modes, however, behave quite differently. Initially, the lower frequency 
mode is a predominantly second rotor-flap mode and the higher-frequency mode is mostly 
a second fuselage bending along with a little second blade-flap mode. Then, above a 
rotor speed of about 30 rad/sec, the two modes exchange behavior so that the first 
becomes fuselage bending and the second becomes blade flapping. In contrast ,. for the 
coupled rotor-fuselage in-air modes (as for the uncoupled modes) , a mode that begins 
as a fuselage-bending mode or as a blade- flapping mode retains its basic character as 
rotor speed is increased, even through a coalescence. 


Coupled Rotor-Fuselage Response 

The vibration response of the coupled rotor-fuselage system in air was investi- 
gated by introducing an arbitrary, uniformly distributed harmonic load along the 
blades (fig. 9) and then determining the response at selected node points on the 
fuselage. It should be noted that the frequency of this loading is not an integer 
multiple of b^^ as it is for a b-b laded rotor rotating at Q. Rather, in order to 
study the dynamic coupling between the rotor and fuselage, all frequencies from 0 to 
120 rad/sec for all rotor speeds from 0 to 60 rad/sec have been considered. The 
absolute magnitude of the acceleration response, in g’s, at the center of the fuselage 
(node 17) owing to a 4,448-N (1,000-lb) harmonic load uniformly distributed along the 
rotor blades is presented as a function of both forcing frequency and rotor speed in 
the 3-D plot of figure 10(a). These same responses are also presented in figure 10(b) 
as a contour map of constant vibration amplitude superimposed on the coupled rotor- 
fuselage eigenvalues and nP excitation frequency lines from figure 5(c). Fig- 
ures 10(a) and 10(b) show that the system has resonant responses at the two lightly 
damped fuselage-bending- mode frequencies (at about 25 and 80 Hz) . Of particular 
interest is the sensitivity of this vibration response to rotor speed near the first 
and second fuselage-bending-mode frequencies. The dynamic coupling between the rotor 
and fuselage modes produces "valleys” in the response topography. The acceleration 
response at the rotor modal frequencies, however, is only significant at low rotor 
speeds, since aerodynamic damping, which increases with rotor speed, suppresses the 
response of the blades and, hence, the fuselage. On the other hand, if aerodynamic 
damping were to be neglected, a considerably different 3-D response plot from that of 
figure 10(a) would be obtained (fig. 11). In this case, the system response at the 
rotor-mode frequencies is not attenuated with increased rotor speed as it is in the 
aerodynamically damped case. The effects of coupled rotor- fusel age dynamics, however, 
are still very much evident. 

The maximum system responses shown in figure 10(a) can be confirmed to coincide 
closely with the coupling of the fuselage- and rotor-mode frequencies by referring to 
the contour-eigenvalue plot in figure 10(b). The minimum response of the second 
fuselage mode occurs at a rotor speed that is roughly equidistant from the rotor 
speeds at which the second and third flap modes coalesce with the second fuselage 
mode. 


If the dynamic coupling between the rotor and fuselage modes were neglected, the 
response would differ greatly from that shown in figure 10(a). Figure 12 shows the 
acceleration response at node 17 for a 4,448-N (1,000-lb) oscillating point-load 
applied at the midpoint of an isolated free-free fuselage. Here the response is seen 
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to vary with forcing frequency o) but, of course, it does not vary with rotor speed. 
In figure 13, the peak responses of this uncoupled fuselage are seen to occur at the 
two uncoupled fuselage symmetric bending-mode frequencies, 25.8 and 83.4 rad /sec. 
Besides these two frequencies at which resonance occurs, there are two other frequen- 
cies worth noting in figure 13: 22.3 and 61.0 rad/sec. At these two frequencies the 

response at node 17 is nearly zero. (If structural damping had not been introduced 
into the problem the response would be exactly zero.) These frequencies are the 
first two natural antiresonant frequencies for node 17 for forcing at node 17. 

Natural antiresonances, which occur in structures much as resonances do (ref. 7), are 
a function of fuselage location, as well as frequency (fig. 14). It is interesting 
to note that the fuselage antiresonant lines in figure 14 represent both the uncoupled 
fuselage forced at its center (node 17) and the coupled rotor-fuselage with forces 
applied to the rotor which is connected to the fuselage at node 17, Although anti- 
resonant theory is not considered directly in the present paper, it is important as 
a possible means of reducing forced vibratory response in helicopters. 

The fuselage nodal responses presented herein represent the absolute magnitude 
of what is a complex response. For completeness, the real and imaginary parts of the 
complex response at nodes 14 and 17 are presented in figures 15(a)-15(d) for the 
coupled rotor-fuselage in air. 

The response plots presented above are for only one particular fuselage node 
point. Since the fuselage is a flexible structure, the response at other fuselage 
node points can be expected to be different. Figures 16(a)-16(f) are included to 
show the variation of the fuselage vibration response from the fuselage tip (node 12 ) 
to the center of the fuselage (node 17) for the coupled rotor-fuselage system sub- 
jected to a 4,448-N (1,000-lb), radially uniform, harmonic load applied to the rotor. 
The response is seen to be strongly dependent on the fuselage nodal location. 

Nodes 12, 14, and 17 are responsive in decreasing magnitude in both fuselage bending 
modes, and nodes 13 and 15 are only responsive in the first and second fuselage- 
bending modes, respectively; node 16 is not very responsive in either mode. These 
modal response results coincide directly with the nodal displacements of the first 
two fuselage mode shapes (figs. 6 (c) and 6 (d)). 


Responses at n per Rev 

It is recognized that in an actual rotorcraft the periodic aerodynamic forces 
generated in forward flight by a rotor with b blades produce fuselage vibratory 
responses only at xibQ (where n = 1,2,3, . • .) since only nbl^ shears are trans- 
mitted from the rotor to the fuselage for a b-bladed rotor. Although response cal- 
culations have been carried out for all combinations of forcing frequency and rotor 
speed, the resonant responses of the system can only occur along the nP = xihQ **per- 
rev” lines on the eigenvalue plots of figures 5(a)-5(c). Thus, for a two-bladed 
rotor only the response along the 2P and 4P lines are meaningful, while for a three- 
or four-b laded rotor only the 3P and 4P lines, respectively, have meaning. The 
higher values of nP are not included here, because the responses diminish rapidly 
as P increases. 

The nP response plot for the coupled rotor-fuselage response at node 17 pre- 
sented previously (fig. 10(b)) is shown in figure 17. In general, the nP response 
peaks occur when the nP lines in figure 10(b) cross a fuselage or blade modal fre- 
quency. In figure 17, the rotor speeds at which the nP lines cross the first and 
second blade-mode frequencies and the first, second, and third fuselage-mode frequen- 
cies are labeled B-^, B 2 , F^ , F 2 , and F 3 , respectively. The largest 2P response 
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peak is seen to occur at a rotor speed of about 40 rad/sec where the 2P excitation 
frequency coincides with the second fuselage-bending-mode frequency (F 2 ) of 
80.7 rad/sec. There are instances, however, where there is no resonant response peak, 
even though an nP line crosses a natural frequency of the system. For example, the 
3P line crosses the second blade-mode frequency (B 2 ) at approximately 18 rad/sec 
(fig. 10(b)), but figure 17 shows no 3P response peak at this rotor speed. In fact, 
the 3P response is very low at this rotor speed. This is due in part to rotor aero- 
dynamic damping but is mostly a result of the fact that the B 2 mode frequency is 
well separated from the Fi and F 2 mode frequencies at this rotor speed (see 
fig. 9(b)). The same is also true where the IP line crosses the first blade-mode 
frequency (Bi). 

Figure 17 also shows how the proximity to a coalescence affects the nP 
response. The 2P and 3P resonant responses at the second fuselage frequency (F 2 ), 
which are near the F 2 mode frequency’s coalescence with the second blade mode, are 
considerably larger than the 4P response at the F 2 mode frequency, which is well 
separated from the coalescences of the F 2 mode with the second and third blade 
modes, respectively. In fact, referring back to the contour plots of figure 10(b), 
the 4P line is seen to intersect the second fuselage-bending mode right in a "val- 
ley” between the two response peaks. The 2P and 3P lines, on the other hand, inter- 
sect the second fuselage mode line high up on the sides of the response "mountain." 

All of the response analyses presented in this paper are for an assumed 6% 
structural damping which was introduced to reduce the magnitude of the resonant peaks 
of the response plots. In order to assess the effect of structural damping on the 
coupled rotor-fuselage response, additional response analyses were carried out with 
2% and 10% structural damping. Figure 18 presents direct comparisons of the coupled 
nP responses at node 17 for 2%, 6% (fig. 17), and 10% structural damping. The fig- 
ure shows the significant effect structural damping has on the resonant response 
peaks that occur at the fuselage natural frequencies. 

In summary, although conventional wisdom says that one should avoid all situa- 
tions in which any system natural frequency is near to or coalesces with an excita- 
tion frequency (because very high resonant responses can occur) , the present results 
show that sometimes these situations may not be a problem at all. Figure 19 presents 
a qualitative representation of regions of low, high, and very high coupled rotor- 
fuselage response at system natural frequencies for node 17. As indicated in the 
figure, an excitation at a point where two system frequencies (i.e., fuselage and 
rotor) coalesce will produce very high vibration response, but if a rotor (or fuse- 
lage) mode frequency lies between two fuselage (or rotor) mode frequencies, then the 
effect of an excitation at that point may be substantially attenuated. Furthermore, 
whereas resonant excitation at fuselage-mode natural frequencies has, in general, a 
large effect on fuselage vibration response, the opposite is seen to be true for 
resonant excitation at blade-mode natural frequencies. Except for low rotor speeds, 
where there is little aerodynamic damping, and near the coalescence of two modes, 
resonant excitation at blade-mode natural frequencies does not lead to large fuselage 
responses. 

Since the coupled rotor- fuselage vibration responses were seen to vary from one 
fuselage node to another (figs. 16 (a)-16 (f ) ) , it is to be expected that the nP 
responses at these nodes would also vary between nodes. Figure 20 presents the nP 
responses at node 14. The resonant nP response peaks for this node are seen to 
occur at the same frequencies as for node 17 (fig. 17), but the amplitudes of these 
responses vary considerably. In addition, the antiresonances, which, as noted earlier, 
vary from one node to another, occur at different frequencies. 
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The variation of the nP responses with fuselage location and rotor speed is 
shown in the 3-D plots of figures 2l(a)-21(d). In each case, the largest response 
is at the ends of the fuselage, although the rotor speeds at which the maximum 
responses occur are different. It should be noted that these 3-D nP plots appear 
somewhat distorted because some of the data points in these figures have been obtained 
by linear interpolation and are plotted using straight lines. Also, these plots are 
of magnitudes only and do not correctly show the node points on the fuselage at which 
the response is zero. 


Approximate Uncoupled Response 

In the preceding discussion, the forced-response behavior of the coupled rotor- 
fuselage system, in hover, was obtained for a uniformly distributed harmonic load 
applied along the rotor blades. In practice, however, coupled rotor-fuselage systems 
have only recently begun to be considered in vibration analyses. Such systems have 
historically been approximated by separate uncoupled systems. Therefore, in order to 
explore the accuracy and influence of this approximation and to investigate the 
effect of dynamic coupling between the fuselage and rotor, response calculations were 
also carried out using the conventional uncoupled approximation for the coupled sys- 
tem. The physical elements of the exact and the approximate calculations are illus- 
trated in figure 22. To clarify the process for the reader, figure 22 traces the 
analyses from applied force to hub shear to fuselage transfer function to vibration 
response. For the uncoupled case, the same uniform 4,448-N (1,000-lb) harmonic load 
previously applied to the rotor of the coupled system was applied to an isolated, 
hub-fixed rotor, and the resulting hub reaction forces (Fp^) were calculated. Then 
in additional separate analyses, these rotor hub reactions were applied to the 
uncoupled, free-free fuselage, both with and without an effective rotor mass included. 
The effective rotor mass in these analyses was taken to be equal to the full rotor 
mass. The uncoupled case with an effective rotor mass is typical of the procedure 
that has been followed by the helicopter industry in an attempt to obtain a better 
approximation of the interaction between the rotor and fuselage (ref. 1). 

In figure 22, the fuselage vibration response at node 17 for the coupled system 
is compared with the response obtained for the uncoupled system, both with and without 
rotor mass. The figure shows that although the qualitative behavior of the responses 
appears to be similar, the magnitude of the approximate, uncoupled results, both with 
and without rotor mass, are considerably larger than the "exact*^ fully coupled 
results. The addition of the rotor mass does, however, result in an improvement in 
both the magnitude of the response and the frequency at which the response occurs. 

This improvement occurs even though the transfer function for the uncoupled fuselage 
without rotor mass is identical to the "exact" transfer function for the coupled sys- 
tem; the approximate transfer function for the uncoupled fuselage with rotor mass is 
somewhat lower in both magnitude and frequency because of the increased mass. 

The cause of the disparity between the coupled and the uncoupled responses is 
apparent when the hub reaction force for the isolated, hub-fixed rotor is compared 
with the hub shear for the coupled rotor-fuselage. Figure 22 shows that in the 
coupled case the presence of the fuselage results in a reduction of the hub shear at 
the forcing frequencies of the fuselage bending modes. For the uncoupled case, how- 
ever, the fuselage modes do not contribute to the hub reaction force. 

A direct comparison of the coupled and uncoupled responses at node 17 is shown 
in the nP response plots of figures 23(a)-23(d). Except for the IP plot where 
the coupled response is the largest, the uncoupled cases consistently overpredict the 
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response of the system. For 2P, 3P, and 4P the uncoupled response without rotor 
mass is seen to be considerably greater than the uncoupled response with rotor mass 
which in turn is considerably greater than the coupled response. 

This underprediction of the IP response and overprediction of the 2P, 3P, 
and 4P responses for the uncoupled cases is shown directly in figures 24(a)-24(d) 
where the ratio of the approximate, uncoupled nP response to the exact, coupled 
nP response is plotted. It should be noted that since figures 24(a)-24(d) represent 
ratios of responses and not actual responses, they are applicable for any node point 
on the fuselage. These figures show, for example, that the peak 4P responses for 
the uncoupled system without rotor mass is almost 6 times greater than the 4P 
response for the coupled system. 

In figures 25(a)-25(d), nP plots of the hub reaction force for the isolated, 
hub-fixed rotor and the hub shear for the coupled rotor-fuselage are presented. 

Except for the IP hub shear, the isolated hub-fixed rotor nP reaction forces are 
seen to be in general considerably greater than the coupled rotor-fuselage nP hub 
shears. As noted above, the influence of the fuselage modes, although clearly evi- 
dent for the coupled case, is absent in the uncoupled case. Finally, these approxi- 
mate and exact nP hub shears are compared directly in figures 26(a)-26(d) where the 
ratio of the isolated hub-fixed rotor reaction force to the coupled rotor-fuselage 
hub shear is plotted. 


CONCLUDING REMARKS 


In this report the dynamic coupling between the rotor and the fuselage of an 
analytical helicopter model in hover was investigated using a simplified beam finite 
element model. Eigenanalyses , as well as forced-response analyses, were carried out 
for both coupled and uncoupled rotor-fuselage configurations. The results for a par- 
ticular beam-model representation of a rotorcraft may be summarized as follows: 

1. The uncoupled rotor-fuselage frequencies 'In vacuo are poor approximations to 
the actual coupled frequencies near rotor-fuselage frequency coalescences. With 
rotor aerodynamic damping included, the uncoupled rotor-fuselage modes are better, 
but still not good, approximations to the coupled system complex modes. 

2. The magnitude of the resonant response at the fuselage-mode frequencies is 
highly dependent on their proximity to the blade modal frequencies. When the rotor 
blade and fuselage frequencies coalesce, vibration is greatly increased (relatively), 
whereas away from such coalescences resonant response at the fuselage-mode frequen- 
cies may be of significantly lower magnitude. 

3. Away from fuselage-mode frequencies, resonant excitation at blade-mode 
natural frequencies has a fairly small and often negligible effect on fuselage vibra- 
tion response because of both dynamics and aerodynamic damping. Thus, it is not 
always necessary to avoid the coalescence of a system natural frequency with an 
excitation frequency. 
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4. For a radially imiform, harmonic blade loading, the approximate vibration 
responses obtained from applying the isolated blade hub shears to a free-free fuselage 
are considerably greater than those obtained when the coupled system is analyzed for 
the same blade loading. 


Ames Research Center 

National Aeronautics and Space Administration 
and 

Aeromechanics Laboratory 

AVRADCOM Research and Technology Laboratories 

Moffett Field, California 94035, August 27, 1982 
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APPENDIX A 


DERIVATION OF ELEMENT MATRICES 


For a beam finite element of length £ (fig* 1), the transverse displacement w 
is related to the modal displacements {q} by 

w = [a]{q} (Al) 


where 


[a] 


"l - 3 ?^ + 25^ 
a - + ^3)^ 

3^2 _ 2^3 
(-5^ + 5^)il . 


(A2) 


and 




{q} 





/ 


(A3) 


For a nonrotating beam element undergoing transverse vibrations, the strain 
energy and kinetic energy are 


= 



dx 


(A4) 

(A5) 


Substituting equation (Al) into the total potential energy 


one obtains 


-d 


J* ^ {q}^[a"]^[a"]{q}dx - j 


j {q}^[a]^[a]{q}dx 


= i {q}^[kE]{q} - J {q}’^[m]{q} 


(A6) 

(A7) 


where the element elastic stiffness matrix [kg] and the element mass matrix [m] are 
given by 
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[kg] = r EI[a"]^[a"]dx 

(A8) 

'o 


fZ T 

[m] = J m[a] [a]dx 

(A9) 


0 


In order to consider beams with nonuniform mass and stiffness properties, the 
present analysis allows the element mass and bending stiffness properties to be repre- 
sented by polynomials of up to third order, that is, 

m = m^ + + mg?^ (AlO) 


El = (EI)o + (EI)i5 + (£1)2?^ + (EDgC^ (All) 

Substituting equations (A2) and (AlO) into equation (A9) and equations (A2) and (All) 
into equation (A8) and carrying out the indicated integrations over the element 
length yields the element mass and elastic stiffness matrices given in equations (Bl) 
and (B2) , respectively, of appendix B. Note that in addition to the distributed mass 
represented by equation (A9) , the present analysis also allows concentrated masses to 
be placed at node points. These point masses are added to the appropriate element 
mass matrix diagonal location. 

For a beam of length R rotating about a fixed axis with angular velocity Q 
there is an increase in the beam’s stiffness as a result of the centrifugal forces 
generated. The element strain energy due to this "centrifugal stiffness" is given by 





T(x) 



(A12) 


where T(x), the tension in the beam element of length Z at a distance r from 
the axis of rotation, is 


T(x) = J dx + Tq (A13) 

•'r 

and where 

R 

dx (A14) 


is the tension due to mass radially outboard of the element of interest. Substituting 
equation (Al) into equation (A12) yields 

f {q}^[k^]{q} (A15) 

where the element centrifugal stiffness matrix [k^] can be calculated from 
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T(x)[a’]^[a']dx 


(A16) 


[kj = 


/ 


r+£ 


Equation (A16) can be evaluated using equations (A2) , (AlO) , and (A13) to obtain the 
centrifugal stiffness matrix given in equation (B3) of appendix B for an element with 
a cubic polynomial mass distribution. 

The aerodynamic loads on the rotor blades of the present coupled rotor-fuselage 
system in hover are obtained from simplified quasi-steady strip theory. The blade 
lift per unit span is given by 


L Poo^c (^^r) (w) 


(A17) 


For an element of the blade the virtual work of this applied aerodynamic load is 


6^ 


1 

= - — Q^aO, \ c(r + x)w6w dx 


(A18) 


Then substituting equation (Al) into equation (A18) yields 


6^ = -{q}*^[c] {6q} 


(A19) 


where 


[c] = ^ c(r + x) [a]^[a]dx 


(A20) 


is the element aerodynamic damping matrix which is presented in equation (B4) of 
appendix B, 

The virtual work of an applied harmonic load p is 


= J p e-^“^6w dx 


£ 


,iu3t i 


(A21) 


Substituting equation (Al) into equation (A21) one obtains 


<5*^ = [p]{<Sq} 


(A22) 


where the element load matrix [p] is given by 


r 


[p] = J p[a]dx e 
0 


io)t 


(A23) 


Equation (A23) is evaluated for a uniform sinusoidal load distribution 
(p = Pq = constant) in equation (B5) of appendix B. 
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APPENDIX B 


ELEMENT MATRICES 


The element mass matrix is 


where 


[m] = 


mil 


sym 

™21 

m 2 2 


®3 1 

”^3 2 

™3 3 

>1 

®4 2 

®4 3 


m 


m 


44 J 


_ 13 3 ^19 ^11 

11 35 + 35 ®i + 630 ”2 s 40 

315 


_ / 11 ^1 ^ 17 1 \, 

"^21 ^210 '*■ 60 '“i 2520 ”^2 315 


-9 ^9 23 

™3i 70 ^ 140 "'i 630 “2 




630 “2 840 “3 


m 


41 


m 


= -(: 
■( 


13 ^1 . 19 . 11 \. 

l^20 “0 70 “1 2520 ™2 2520 ^3)^ 


2520 3^ 

1 ^1 ^1 ^ 1 \, 
’^n 3 QO + 630 "*2 1260 ”^3/ 


^22 ^105 “0 180 ‘“1 

1 


-/13 ^1 ,..5 _^2 \„ 

“32 ^420 “0 160 504 “2 + 335 

“42" "(420 ®o ■*■ 2M ™i 5^ “2 840 


m 0 s — 


13 


29 


23 


‘33 35 ™o 7 ®i '*’ 326 ™2 "*■ 320 


m 


4 3 


_/ll ^1 ^13 ^7 \. 

\210 "^0 28 “1 504 ”*2 ■'■ 360 


“44 (305 “0 368 ™i 252 “2 ■'■ 360 “3)^ 


(Bl) 
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The element elastic stiffness matrix is 


where 




^11 


^2 1 ‘''2 2 


^31 ^32 


L^4l 


'•4 2 


symmetric 


"33 


^3 


"11 


^21 


*"3 1 


"41 


"2 2 


"32 


"4 2 


"33 


[12(EI)o + 6 (EI)^ + 4.8(EI>2 + 4 . 2 (El) 3 ] /«,^ 
[ 6 (EI )3 + 2(EI)^ + 1.4(EI)2 + 1.2(EI)3]/J12 
-1^11 

[6(EI)o + 4(EI)i + 3.4(EI)2 + 3(EI)3]/£2 
[4(EI)„ + (EI)^ + (8/15)(EI)2 + 4(EI)3]/JI 
-^21 

[2(EI)o + (EI)^ + (13/15) (ED 2 + 8 (EI) 3 ]/«, 

k. 


"11 


^43 ^41 


= [4(ED^ + 3(EI)^ + (38/15)(EI)^ + 2.2(EI)J/Jl 


The element centrifugal stiffness matrix is 


(B2) 


(t "^0 + -J "12 + 3 "*3) + ra^^mg + ^ ^ m2 + -| m3^|[k^^] 

- m^rQ^tkc^] - j (m^Z + mj^r)J5^ [k^. ^ ^ (m^Jl + m2r)S^^[kc ] 


- (m2«. + m3r)S^^[k^2 ] - (m3«,)J2^ [k^, ] 


(B3) 
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where 



"6/5 


symmetric 



i/10 

11^115 



[kcj = 

-6/5 

-i/10 

6/5 



i/10 

-i^/30 

-i/10 2i^/15j 


"3/5 


- 





symmetric 



i/10 

i\/30 



] = 




2 

-3/5 

-i/10 

3/5 



0 

-iVeo 

0 i^/10_ 



"12/35 







symmetric 


i/14 

2i^/105 



[l^cj = 






-12/35 

-i/14 

12/35 



-£/35 - 1^/10 £/35 3£^/3Sj 


[k, ] 


"3/14 


i/20 

lli^/840 

-3/14 

-i/20 

-i/28 

-lli^/840 


symmetric 

3/14 

£/28 13£^/168 


[k, ] 


"l/7 

i/28 

-1/7 

-i/28 


symmetric 

i^/105 

-i/28 1/7 

-i^/84 i/28 i^/14. 


[kc ] 


" 1/10 
lli/420 

i^/140 

S3mimetric 

-1/10 

-lli/420 

1/10 

_-i/30 

-3i^/280 

i/30 i^/15. 
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To = 


in n 


m 

[r 2 - (r + <0^] +_i [r3 _ (r + £)3] 


m 


+ ^ ^ [R** - (r + £)"] [r 5 - (r + £)^] 


m 


3 1 


^ £2 


5 


The element aerodynamic damping 'matrix is 


11 L + A. 

35 £ 35 


[c] = j p^ac^A 


( 1 Al + a_v + -i^£^ 

\210 £ 60 ) \105 £ 280^ 


symmetric 


9 r 9 


( 11 - 


+ 


70 £ 140 \420 £ 60 


my 


11 L + 1 

35 £ 7 


[_ (420 £ ^ 70 )^ (l40 £ 28o)^^ ( 2^0 £ 2^)^ (l05 £ 168)*' J 


(B4) 


The element load matrix [p] for a uniform harmonic load is given by 

/ 


T 

[P] = P. 


1 


1/2 

£^/12 

£/2 

-£^/12 


i03t 


(B5) 
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APPENDIX C 


EQUATIONS OF MOTION 


The equations of motion for an element of a rotating beam can be derived from 
Hamilton’s equation, 

t2 

- J 6^dt = 0 (Cl) 

1 

Substituting equations (A7), (A15), (A19), and (A22) of appendix A into equation (Cl) 
and carrying out the indicated variations yield the element equations of motion, 

[m]{q> + [c]{q> + ([kg] + [k^]){q> = {p} (C2) 

For a system made up of more than one finite element, the element mass, damping, 
and stiffness matrices [m] , [c] , and ([kg] + [k^]) and element load vector {p} can be 
assembled into the corresponding global matrices [M] , [C] , and ([Kg] + [K^]), and 
global vector {P}. Then equation (C2) may be rewritten as 

[M]{q} + [C]{q> + ([Kg] + [Kc]){q} = {P> (C3) 



This is the basic matrix equation of motion considered in the present analysis. 

For the case of free harmonic vibration (in vacuo) , the damping matrix [C] and 
the load vector {P} are set to zero in equation (C3) and we obtain the eigenvalue 
problem. 


rX J[M]{q} = ([Kg] + [Ke]){q} (C4) 

which has the real eigenvalues X == con and the matrix of real eigenvectors [U] . 

In carrying out the complex eigenvalue analysis of equation (C3) , it is useful 
to uncouple the equations of motion as outlined by Meirovitch (ref. 8). A coordinate 
transformation is first performed on equation (C3) by replacing 

{q} = [UHi) (C5) 

where [U] is the matrix of orthonormal modes associated with the real eigenvalue 
analysis of equation (C4) . Then equation (C3) is premultiplied by [U]^, yielding 

[U]^[M][U]{q} + [U]^[C][U]{t> + [U]^([Kg] + [Kj,])[U]{q} = [U]^{p} (C6) 

Since the normal modes [U] from equation (C4) are such that 

[U]’’^[M][U] = [I] (C7) 

and 

[U]^([Ke 1 + [Kj,])[U] = ro)2 J (C8) 


20 



where [I] is the identity matrix, equation (C6) becomes 

[I]{q} + [C]{q} + r^Wn'J{q} = {P} 

where 

[C] = [U]^[C][U] 

and 

{P} = [U]'^{P} 

The complex eigenvalues of equation (C9) can then be determined from 


/i 

q 


[ 0 ] 




[I] 

-[C] 


|q I 

|t| 


Alternatively, for harmonic excitations. 


{P} - 

and 

{q} = 


so that equation (C9) can be rewritten as 
where 

{Pq} = tu]'^{p„} 


(C9) 

(CIO) 

(Cll) 

(C12) 

(C13) 

(C14) 

(CIS) 


The responses obtained from equation (C15) are the modal responses as a function 
of forcing frequency oo for a given applied load vector {Pq}. The corresponding 
nodal responses {qo) can then be obtained from equation (C5). 

In the frequency-response analysis, structural damping can be readily included 
by introducing the complex stiffness k(l + ig) in place of k where g > 0 is the 
structural damping factor (ref. 9). Then equation (C15) becomes 

+ ig) - + i‘^[C]){qo> = {Pq} (C16) 

In the present finite element analysis, the total number of degrees of freedom 
is equal to twice the number of node points minus the number of rigid body constraints. 
In general, however, only a relatively small number of these modes are of interest. 

In this case, the formulation described above can be used to obtain approximate, but 
still very accurate results if a reduced set of only the lower frequency modes neces- 
sary to give sufficiently converged results is used in equation (C6) . This approxi- 
mation has the advantage that instead of solving the complete system of equations 
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presented in equation (C9) one can achieve nearly as accurate results for the lower- 
frequency modes of interest with a reduced system of equations. All of the results 
presented in this report, however, were obtained from analyses which include all of 
the collective rotor modes and symmetric fuselage- bending modes. 
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Figure 1.- Be am- element geometry. 
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Figure 2.- Rotorcraft model. 



El, N-m^ X 10® ^ MASS/LENGTH, kg/m 


POINT MASS (NODE 6) = 233.5 kg (16.0 slug) 



-8-4 0 4 

X, m 


■igure 3.- Rotor and fuselage mass distributions. 
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Figure 4.- Rotor and fuselage stiffness distributions 
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(a) In vaouo. 

Figure 5.- Coupled and xincoupled rotor-fuselage eigenvalues 










(a) First rotor-flap mode. 

Figure 6.- Uncoupled and coupled rotor-fuselage mode shapes. 
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(d) Second symmetric fuselage-bending mode 


Figure 6.- Concluded. 
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(a) Second rotor-flap mode, 9 . = 24-36 rad/sec. 

Figure 7,- Coupled rotor-fuselage mode shapes near coalescences. 
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(b) Second symmetric fuselage-bent 


Figure 7.- 


mode, Q = 24-36 rad/sec 


ntinued. 







mode, = 22-28 rad/sec. 
- Continued. 





(d) First symmetric fuselage-bending mode, Q = 22-28 rad/sec 

Figure 7.- Concluded. 






Figure 8.- Effect of dynamic coupling and aerodynamic damping on the 
coupled rotor-fuselage modal frequencies near a coalescence. 
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F = sin cjt 



Figure 9.- Excitation for coupled rotor-fuselage model. 
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Figure 10.- Coupled rotor-fut 
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Figure 12.- Uncoupled fuselage response at node 17 to an oscillatory point-load at 
rotor hub location; acceleration versus forcing frequency and rotor speed. 
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(d) Node 15 response. 
Figure 16,- Continued. 






(f) Node 17 response. 
Figure 16.- Concluded. 
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Figure 19.- Qualitative representation of coupled rotor-fuselage response at 

system natural frequencies at node 17. 
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Figure 24.- Ratio of approximate, uncoupled nP response to exact coupled nP response 
(a) IP response; (b) 2P response; (c) 3P response; (d) 4P response. 
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